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Abstract 

Neutrino mass spectrum is reanalyzed in supersymmetric models with 
explicit trilinear R violation. Models in this category are argued to provide 
simultaneous solution to the solar and atmospheric neutrino anomalies. It 
is shown specifically that large mixing and hierarchical masses needed for 
the vacuum solution of neutrino anomalies arise naturally in these models 
without requiring any additional symmetries or hierarchies among the trilinear 
couplings. 
1. Introduction: 

The hypothesis of neutrino oscillations has gained acceptance after careful observation of 
the atmospheric muon neutrino deficit at the Superkamioka If neutrinos do oscillate 
then both the solar and the atmospheric neutrino deficits can be simultaneously understood 
in terms of vacuum oscillations among the three known neutrinos. This however requires 
the presence of two large mixing angles among three neutrino states. Many different 
theoretical models P have been proposed in this context. Supersymmetry (SUSY) provides 
a framework where both the largeness of mixing and hierarchy in masses can be naturally 
understood. 

Supersymmetric extension of the standard model contains the following lepton number 
violating terms: 

= Kjk^iQjDl + XijkLiLjEl + eiLiH2 . (1) 

These naturally lead to neutrino masses ||^. The neutrino spectrum in this model has 
been extensively studied in the literature in recent times. It has been shown 



that bilinear SUSY violating interactions provide very economical framework which can 
simultaneously accommodate hierarchical masses and large mixing. In contrast, the neutrino 
spectrum implied by the trilinear interactions A, A' may appear arbitrary a priori due to very 
large number of such couplings. It was emphasized by Drees et al. [0 that this is not the 
case and neutrino spectrum could be quite predictive even in models with trilinear couplings. 
This can be understood from eq.(|l|) which in the absence of bilinear terms is invariant under 
a global U(l) symmetry with U(l) charges 1, -1, -2 for the fields L, D'^, E'^ respectively. This 
symmetry would thus prevent generation of neutrino masses if it was not broken by the 
down quark and charged lepton Yukawa couplings and /if respectively. This means that 
the neutrino masses generated by the trilinear interactions in eq.(|ID are always accompanied 
by the above Yukawa couplings and hierarchy in the latter gets decoded into the neutrino 
masses if all the trilinear couplings are assumed to be similar in magnitude. In the limit 
of keeping only the b quark Yukawa coupling, one combination of neutrino fields namely 
A^33t'i, obtains a mass and this combination would contain large mixing of all the states if 
A^33 are comparable for different i. The other mass would arise when the strange quark 
Yukawa coupling is turned on and one thus naturally gets [0], 



(2) 
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reproducing the hierarchy needed to understand the solar and atmospheric neutrino anoma- 
lies simultaneously. The natural expectation of this scenario is large mixing among all three 



neutrinos and this is not favored by the more likely small angle MSW ||TlI] solution for the 
solar neutrino deficit. This led to imposition of ad- hoc discrete symmetries in to prevent 
unwanted trilinear couplings reducing the attractiveness of the scenario. 

It is clear from the forgoing discussion that more natural possibility with the trilinear 
couplings of similar magnitudes is to have large mixing among all the neutrinos. This 
however then favors the vacuum solution to the solar neutrino problem in which case the 
hierarchy among neutrino masses is required to be stronger than displayed in eq.(0). A 
careful analysis of the neutrino spectrum reveals that under the standard assumptions, the 
neutrino mass hierarchy resulting in models with only trilinear R violating couplings at 
a high scale is indeed stronger than the one in eq.(0). It can be strong enough to get 



the vacuum solution for the solar neutrino problem. This coupled with large mixing among 
neutrinos alleviates any need to postulate discrete symmetry as in ^ and makes the trilinear 
lepton number violation an attractive means to understand neutrino anomalies. 

The key feature leading to a different conclusion compared to [0 is the observation 
that the presence of trilinear interactions in the original superpotential at a high scale, 



induces [|12[ terms linear in the sneutrino fields in the effective potential at the weak scale. 
These sneutrino fields then obtain vacuum expectation value (vev) and cause neutrino- 
neutralino mixing. Neutrino mass generated through this mixing dominates over the loop 
mass considered in and in other works [^|13|- This alters the neutrino mass hierarchy 
compared to eq.(|D. We discuss these issues quantitatively in the following. 
2. Sneutrino Vevs and Neutrino Masses 

For definiteness, we shall concentrate on the trilinear interactions containing A' couplings 
and comment on the inclusion of the A couplings latter on. The presence of non-zero A^^-^j is 
known to induce two separate contributions to the neutrino masses and we discuss them in 
turn. 

A. Tree level mass 



We adopt the conventional supergravity framework according to which the structure 
of the superpotential dictates the structure of the soft SUSY breaking terms. Thus, with 
only trilinear L- violating interactions, the soft terms do not contain bilinear terms at a high 



scale. They are nevertheless generated at the weak scale |jT2[ and should be retained in the 
scalar potential at this scale: 



B,H',Hl - B^J^.Hl + h.c] + ^-{gl + gl){\ p - | + (3) 

Where, we have retained only neutral fields and used standard notation with B^. and ^tiHi 
representing the bilinear lepton number violating |T5| soft terms. The weak scale value of the 
soft parameters is determined by the following |]T2[ renormalization group equations (RGE): 



dt 

d'^^Hi 2 / 3 , / 3 

— m 



Hkk I 1 



and the standard RGE |TJ] for the parameters on the RHS. Since we allow only trilinear 
interactions in VT^, rn^v-n^ = B^- = at high scale. As seen from the above equations, the 
presence of non-zero X[f.i^ however generate non-zero values for ttlI,^^ and B^-. It is then 
convenient to parameterize them as, 

Bei = ^ipphp Kip, 

Here, p is summed over generations. The parameters k and k' represent the running of the 
parameters present in the RGB's from the GUT scale to the weak scale. 
The above soft potential would now give rise to sneutrino vevs, 

< J^i > = — n 1 — — ■ 6 

m% + \m\ cos2p ^ ' 

The sneutrino vevs so generated will now mix the neutrinos with the neutralinos thus giving 



rise to a tree level neutrino mass matrix [0 : 



2{-cfiM2 + 2M^cpSp{c + tanei)) ' 



B. Loop Level Mass 

The trilinear couplings in the superpotential would also give rise to a loop induced neutrino 
mass with the down squark and antisquark pairs being exchanged in the loops along with 
their ordinary partners This mass can be written as, 

3 

= Y^Xifc^j-fci ^1 sincpi cos(l)i . (8) 

In the above, sirKpi cos(pi determines the mixing of the squark-antisquark pairs and and 
M|; represent the eigenvalues of the standard 2x2 mass matrix of the down squark system 
I^ . The indices / and k are summed over. The Vi stands for the vev of the Higgs field H^. 



The mixing sin(j)i cos(pi is proportional to /iP and thus one can write the loop mass as, 

M\j = Kik >^'jki hk hf mioop . (9) 



Explicitly, 



3^1 sincpi coscpi M| 3 vf 1 
167r^ hi M(i IQtt^ Msusy 



with Msusy ~ 1 TeV referring to the typical scale of SUSY breaking. Note that mioop 
defined above is independent of the R violating couplings and is solely determined by the 
parameters of the minimal supersymmetric standard model (MSSM). 

As evident from eqs.(0,^, the tree as well as loop induced masses have very similar 
structure. Both involve down-quark Yukawa couplings for the reason explained in the intro- 
duction. The tree level contribution involves diagonal couplings A^^;. while the A^' contains 
off-diagonal A-^^ as well [ffSl. If all A-;,; are assumed to be of similar magnitude then the tree 



ikl 

level mass is seen to dominate over the loop mass as we will discuss in the next section. 
3. Neutrino Masses and Mixing 

We now make a simplifying approximation which allows us to discuss neutrino masses and 
mixing analytically. It is seen from the RG eqs.(^) that the parameters Kjfc, k^^, defined in 
eq. d^) are independent of generation structure in the limit in which generation dependence 
of the scalar masses m|., mg. and soft parameters and Af^^j^ is neglected. Since we are 
assuming the universal boundary conditions, this is true in the leading order in which the 
dependence of the parameters multiplying X'lkjji'j^ in eq.® is neglected. dependence 
in these parameters generated through the gauge couplings will also be fiavor blind though 
Yukawa couplings will lead to some generation dependent corrections. But their impact 
on the the conclusions based on the analytic approximation below is not expected to be 
significant. The neglect of the generation dependence of KjA;, i^'ik allows us to rewrite eq.(^ 
as, 

M% = nioaiaj , (11) 

where, 

a^ ^ Kkk (12) 

rriQ is now completely determined by the standard MSSM parameters and the dependence of 
the R- violating parameters gets factored out as in eq.(j^). mo can be determined by solving 
the RGE (§). Roughly, mo is given by. 



3 \2 

Let us rewrite the loop induced mass matrix as, 



Adli — rriioop ^iik Kki 



In 



M 



xV 



Ml 



mioop CLi aj + mioop Aij + , ) , 



where. 



A — \' \' _i_ X' \' _ X' X' — X' X' 

^ij ~ ^i23^j32 ' ^i32^j23 '^i22'^j33 ^i33^j22 



(13) 



(14) 



(15) 



Neglecting 0(/if ) corrections to the loop induced mass matrix, the total mass matrix is 
given by. 



^ (mo + mioop) cbittj + miocyp A. 



^3 



M'ij + mioop /if ^ 



(16) 



The matrix M' has a special structure. It has only one eigenvalue. It can be easily diago- 
nalised using a unitary transformation. 



U^M'U = diag (0,0,7713) 



(17) 



where. 



ma a; (mo + mioop) («! + al + aj) 



(mo + mioop) ( ^^333 + -^233 + ^ISd) ^3 



'2 



/2 \ uD^ 



The matrix U is determined as. 



( \ 

C2 S2C3 S2'S3 
-S2 C2C3 C2S3 



y -S3 C3 y 



with. 



S2 



oi 



S3 



(18) 



(19) 



(20) 



The total mass matrix is now given by, 



U^M^U ^ mg diag (0, 0, 1) + niioop hf A' 



( 



eA[^ eA[2 eA[^ 



^^22 ^^23 



eA[^ 6^23 



(21) 



where, 



A' = f/^A U 



(22) 



and 



^ A' ~ ''^loop r,D r^D Al ^ ^loop ^2 

mo fts 



The last equality follows under the assumption that A^^-^ are similar in magnitude and 
Choosing, 



U' 



I \ 

c\ s\ 

-S\ Ci 

1 



(23) 



with, si, ci defined by. 



tan 2^1- ^^^'^ 



A' — A' ' 

^22 ^11 



(24) 



we have. 



U'^U^M'^UU' = ms 









e52 



y Ci€A[^ — Sie/1.23 — CieA: 



23 



sieA[s - cieA'23 
1 



ma diag{e6i, 6^2, 1) • 



(25) 



The off-diagonal elements will generate additional mixing in the model. But, as eA' <^ 1, 
we can neglect these off-diagonal elements. The eigenvalues in this approximation are given 
as. 



m 



1^1 



(26) 



where, 



5i = {cl - 2ciSiA[2 + sjA 



2 /I/ \ 

22) > 



{sf A[^ + 2ciSiA[2 + cjA' 



22 J 



(27) 



Note that both 6i and 62 are generically of 0(A ) when all A^^^ are assumed to be similar 
in magnitude. As a consequence, the neutrino masses follow the hierarchy. 



With, 



'^s mi 
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m^3 rUb niQ V^i-^i33, 



(28) 



The last factor in the above is of 0(1) and the remaining part is controlled completely by 
the standard parameters of the MSSM. 

Eq. may be regarded as a generic prediction of the model. It is seen from eqs. (|TD|,|T^ 
that typically. 



mi 



oop 



TT 



mo 
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10 



-3 



(29) 



Thus the neutrino mass ratio in eq.(|28|) is suppressed considerably compared to eq.(0) ob- 
tained when sneutrino vev contribution is completely neglected. The exact value of this 
suppression factor is dependent on MSSM parameters and we will calculate it in the next 
section. 

The mixing among neutrinos is governed by. 



K = U U' 

/ \ 

C1C2 - S1S2C3 S1C2 + C1S2C3 S2S3 

-S2C1 - S1C2C3 -S1S2 + C1C2C3 C2S3 



S1S3 



-S3C1 



C3 



(30) 



The angles are determined by the ratios of the trihnear couphngs and hence can be natu- 
rally large. Thus, as in supersymmetric model with purely bilinear R violation one gets 
hierarchical masses and large mixing without fine tuning in any parameters. 
4. Neutrino Anomalies 

We now discuss the phenomenological implications of neutrino masses, eq.(PB|) and mixing, 



eq. (|30D . Due to hierarchy in masses, one could simultaneously solve the solar and atmo- 
spheric u problems provided, rriui ~ ''^^2 ~ 10"^ eV and m,i^.j ~ 10~^ eV. 

In order to determine these masses exactly, we have numerically integrated eqs.@ along 
with similar equations for the parameters appearing in them. We have imposed the standard 
universal boundary condition and required radiative breaking of the SU{2) x U{1) symmetry. 
Solution of the RGE determines both rriioop (eq.®) and m.Q (eq.(pTD). We display these in 
figs, (la, lb) as a function of /i for tan(3 = 2.1, M2 = 200 and 400 GeV respectively. The 
ratio is quite sensitive to the sign of For /i > 0, this ratio is rather small, typically, 
~ 10^^ — 10^'^, while it can be much larger for /i < 0. There exists a region with negative /i 
in which > 1. In this region, two contributions to the sneutrino vev in eq. (|^) cancel 
and mo gets suppressed. Barring this region, the ^^^^ is seen to be around ~ 10^^ — 10^^ 



for negative fj, leading to 



m,, _ nismi^op . , 3 (iQ-^ _ IQ-^) (31) 



rsj r\j 



For m^yg ~ 10~^ — 10~^ eV, one thus obtains m^^ ~ ''^ui ~ 2 (10~^ — 10~^) eV which is in the 
right range required to solve the solar neutrino problem through vacuum oscillations. The 
typical value of ttiq ~ GeV found in Fig. (lb) implies through eq. (p!8D , 

A' ~ 10-^ 

Thus one needs to choose all A^^j^. of this order. Once this is done, one automatically obtains 
solar neutrino scale for some range in the MSSM parameters. 

While, hierarchy needed for the vacuum solution follows more naturally, one could also 
obtain scales relevant to the MSW conversion. This happens for very specific region of 
parameters with negative fi in which two contributions to sneutrino vev, eq. (P), cancel. As 
already mentioned, can be 1 in this region. One then recovers the result of , namely. 



eq.(0) which allows MSW solution for the solar neutrino problem. The reference P] which 
used hierarchical A^^^ also concentrated on this region in order to obtain the MSW solution. 

We showed in Fig. 1 neutrino mass ratio for specific value of M2 and tan/?. Qualitatively 
similar results follow for other values of these parameters. We have displayed in Table. 1 
values for the MSSM parameters and what they imply for We have shown illustrative 
values of the parameters which lead to the vacuum as well as MSW solution. The latter 
arise only for limited parameter range corresponding to cancellations in eq.@. The former 
is a more generic possibility which arise for larger region with both positive and negative 
values of /i. The MSW solution in the present context will have to be restricted to the large 
angle solution if one does not want to impose any discrete symmetries or fine tune A's. 

The constraints on mixing matrix K, eq.(0), implied by the experimental results are 
also easy to satisfy keeping all the A^^^ similar in magnitude. Hierarchy in masses, rriy^, m^^ 
lead to the following survival probabilities for the solar and atmospheric neutrinos after 
undergoing vacuum oscillations: 

= 1 - 4 K',, Kl, sin^ - 2 A- (1 - A- ) (32) 

Where, Am^^ = m^^ — rrv^. and Am|;^ ~ Amgg. 

These survival probabilities assume the standard two generation form when K^^i = and 
one could utilize existing constraints on mixing angles. In practice, the K^^ may not be zero 
but is constrained to be small from the non-observation fl^ of z/g oscillations at CHOOZ. 
To the extent it is small, one could use the two generation constraints for the solar and 
atmospheric analysis. The combined constraints which are needed [0,^ to be satisfied are: 

0.6 < AKl^ = si, sin2 2^2 + 4 sin^ 2^3 < 1. 

Ke3 < 0.18 (33) 

0.8 < 4 K^^, Kl^ = A{ciC2 - sis^c^f {s^c^ + s^cic^f < 1. 

It is possible to satisfy all these constraints by choosing for example, 

C3 = S3 = si = ci = ; S2 = 0.28 



The relative smallness of S2 required here does not imply significant fine tuning and can be 
easily obtained, e.g. by choosing, 

■^133 1 
-^233 3 

We have so far concentrated on the A^^^ couplings alone. The analogous discussion can 
be carried out for Xij^ couplings appearing in the eq.(0). Here also, the tree level contribu- 
tion to neutrino masses will dominate over the loop contribution although the structure of 
mixing matrix will differ slightly due to the anti-symmetry of the couplings Xijk in indices i 
and j. 

5. Discussion : 

We have discussed in detail the structure of neutrino masses and mixing in MSSM in the pres- 
ence of trilinear R- violating couplings, specifically X'^jf.. Noteworthy feature of the present 
analysis is that it is possible to obtain the required neutrino mass pattern under fairly gen- 
eral assumption of all the X[jf^ being of equal magnitudes. This is to be contrasted with 
the recent analysis [0-^] which had to make very specific choice of the trilinear couplings 
in order to reproduce neutrino mass pattern. It is quite interesting that hierarchy among 
neutrino masses is controlled by few parameters in MSSM and is largely independent of the 
trilinear R violating couplings. Thus one could understand the required neutrino mass ratio 
without being specific about the exact values of large number of the trilinear couplings. This 
'model-independence' is an attractive feature of the scenario discussed here. 

The key difference of the present work compared to many of the other works is proper 
inclusion of the sneutrino vev contribution. While we had to resort to specific case of the 
minimal supergravity model for calculational purpose, the sneutrino vev contribution would 
arise in any other scheme with A^^^ 7^ at a high scale such as Mqut- Such contribution 
thus cannot be neglected a priori. On the contrary, the inclusion of this contribution makes 
the model more interesting and fairly predictive in spite of the presence of large number of 
unknown couplings. 
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Tablel: The values of mo, mioop sind ratio of the eigenvalues, for various values of the 
standard MSSM parameters m, M2 and ji for tan (3 — 2.1 , A—0. 
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Figure la. The absolute values of tree level contribution, mo, the loop level contribution, 
mioop and their ratio ^^^^ are plotted with respect to ji (positive) for M2 — 200 GeV, A—0 
and tanP — 2.1 . The mo and mioop are defined in the text. 
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Figure lb. Same as in Fig. (la) but for — 400 GeV and (negative). 
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